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Is There a Planet in the PSR 1620-26 Triple System? 

Frederic A. Rasio 1 
Institute for Advanced Study, Olden Lane, Princeton, NJ 08540 

ABSTRACT 

The unusually large eccentricity (ei = 0.025) of the low-mass binary millisecond 
pulsar PSR B1620-26 can be explained naturally as arising from the secular perturbation 
of a second, more distant companion. Such a triple configuration has been proposed 
recently as the most likely cause of the anomalous second period derivative of the 
pulsar. The current timing data are consistent with a second companion mass rri2 as 
low as ~ 1O _3 M0, i.e., comparable to that of Jupiter. However, if the eccentricity is 
indeed produced by secular perturbations, then the second companion must be another 
star, most likely of mass rri2 <J IMq and in a very eccentric {e.2 i; 0.5) orbit of period 
P2 ~ 10 2 -10 3 yr. A second companion of planetary mass cannot induce the observed 
eccentricity. Independent of the mass of the second companion, small changes in the 
binary pulsar's orbit should become detectable with just a few additional years of timing 
data. This detection would provide direct confirmation of the triple nature of the system, 
and an accurate measurement of the effects would place important new constraints on 
the orbital parameters. 

Subject headings: celestial mechanics, stellar dynamics — planetary systems — pulsars: 
general — pulsars: individual (PSR B1620-26) 

1. Introduction 

The millisecond pulsar PSR B1620-26 in the globular cluster M4 has a low-mass companion 
(mass mi ~ O.3M for a pulsar mass m p = 1.35 M & ) in a nearly circular orbit of period Pi = 
0.524 yr (Lyne et al. 1988; McKenna & Lyne 1988). The eccentricity e x = 0.0253, although small, is 
several orders of magnitude larger than observed in most other low-mass binary millisecond pulsars 
(Thorsett, Arzoumanian, & Taylor 1993, hereafter TAT). In addition, the timing data indicate 
the presence of a very large second derivative of the pulse period, P = —2.3 X 10 _27 s _1 (Backer 
1993). This is many orders of magnitude larger than would be expected from a varying acceleration 
caused by either the mean gravitational field of M4 or a nearby passing star. Instead, the most 
likely explanation is that the pulsar has a second, more distant orbital companion (Backer 1993; 
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Backer, Foster, & Sallmen 1993; TAT). Such hierarchical triple systems are expected to be produced 
quite easily in globular clusters through dynamical interactions between binaries (Mikkola 1984; 
Hut 1992). 

The present timing data are consistent with a second companion mass anywhere in the range 
10~ 3 M & <J rri2 <J l-^© with a corresponding orbital period 10 yr <J P2 <J 10 3 yr. This opens 
the exciting possibility that the second companion could be a Jupiter-like planet (Backer 1993; 
Sigurdsson 1993; TAT), although it could also be another star (main sequence, white dwarf, or even 
another neutron star). Planets have been detected in orbit around at least one other millisecond 
pulsar, PSR B1257+12 (Wolszczan & Frail 1992). 

If PSR B1620-26 is indeed in a triple system, then small perturbations of the (inner) binary 
pulsar's orbit should be induced by the presence of the second (more distant) companion. In §2 
below, we show that the unusually large eccentricity of the binary pulsar can be explained naturally 
as arising from these perturbations, but only if the second companion is a star and not a planet. 
In either case, small changes in the eccentricity and longitude of pericenter of the binary pulsar 
should become detectable with a few additional years of timing data (§3). Such a detection would 
provide independent confirmation that the first pulsar in a hierarchical triple system has indeed 
been discovered. The situation is similar to that of PSR B1257+12, where the recent detection of 
orbital perturbation effects has confirmed the presence of (at least) two planets in orbit around the 
pulsar (Rasio et al. 1992; Wolszczan 1994). 

2. Secular Perturbations 

Assume that the binary pulsar had a very low eccentricity ei(0) <C = 0.025 at the 

time it acquired its second companion. Standard secular perturbation theory of orbital mechanics 
(e.g., Brouwer & Clemence 1961) can be used to calculate the subsequent evolution of the system 
over a time t p >• P2 > P\. The particular solution for the eccentricity ei(t) of the inner orbit 
corresponding to an initial condition with ei(0) = takes the form 

ei (f) = F(q u q 2 , a) e 2 (0) [1 - cos(^)] 1 / 2 (1) 

Here J 7 is a dimensionless function of qj = mj/m p and a = 0,1/0,2 < 1. The angular frequency 
g = <7i — <72j where g\ and 52 are the eigenvalues of the secular perturbation matrix. An expression 
for T can be written explicitly in terms of Laplace coefficients b^ 2 (a) (Brouwer & Clemence 1961; 
a useful summary of all relevant expressions is given by Dermott & Nicholson 1986). For a<l 
and q^/qi >■ a 1 / 2 we find T ~ a (independent of q?) and 2ir / g ~ Piq^ 1 a~ 3 , whereas for a <C 1 
and 92/91 *C a 1 / 2 , T ~ q^a 1 ! 2 and 2ir / g ~ Pia -7 / 2 . Other treatments of secular evolution in triple 
systems are given by Mazeh & Shaham (1979), and Bailyn (1987). 

The analytic solution (1) is valid only for small eccentricities and small relative inclination. 
If the system was formed through a dynamical interaction, the eccentricity and inclination of 
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the outer orbit are likely to be quite large. Numerical integrations of the three-body problem 
must then be used. We have constructed a large number of numerical solutions covering the full 
range of parameter values allowed by the current timing data (Rasio 1994). We find that the 
form of the solution remains generally similar to that of equation (1), with e\{t) returning to zero 
periodically every ~ 10 3 -10 6 yr. However, both the period and amplitude of the solutions can 
deviate significantly from their analytic values when e 2 or the inclination is large. 

Our most important result is illustrated in figures 1 and 2. When the mass of the second 
companion is ^ 0.1 Mq, the maximum eccentricity (Aei) max induced in the binary pulsar's orbit 
remains always smaller than the value e\ = 0.025 observed today. For a typical planetary case, 
with m 2 = 80 M ffi and P 2 = 10 yr (cf. TAT), we get (Aei) max ^ 0.001 for all e 2 . In contrast, for a 
typical stellar case, with m 2 = O.8M and P 2 = 120 yr (TAT), we find that (Aei) max > 0.025 can 
be obtained easily for e 2 > 0.4. We conclude that, if the eccentricity of the binary pulsar's orbit 
has been induced by secular perturbations, then the second companion must be a star rather than a 
planet. 

Note that, according to equation (1), a larger amplitude (Aei) max can always be produced 
by invoking a larger eccentricity e 2 for the outer orbit. However, dynamical stability (e.g., Black 
1982) requires that (1 — e 2 )a 2 /ai ^ 3-4, placing a direct limit on how large e 2 can be for a given a 
(cf. Fig. 2). For example, we note that the solution with m 2 = 0.8 M@, P 2 = 120 yr, and e 2 = 0.9 
mentioned by TAT is dynamically unstable. 

A number of alternative explanations for the anomalous e\ are possible, depending on the triple 
formation process. If the second companion is a star, then the most natural formation process is 
through an interaction between the (pre-existing) binary pulsar and a much wider primordial binary 
containing two main-sequence stars. The typical outcome of such an interaction is a triple system 
containing one of the main-sequence stars in a wide orbit around the binary pulsar. The orbital 
eccentricity of the binary pulsar could have been perturbed during this interaction. However, this 
would require a close passage (within r <J 5ai) by one of the main-sequence stars, implying a high 
probability of disrupting the binary pulsar completely. A later close encounter with a passing star 
in the cluster is even less likely, since it has a high probability of also ejecting the second companion. 

Sigurdsson (1993) describes a formation scenario where a pre-existing neutron-star- white-dwarf 
binary has an interaction with a turn-off main-sequence star with a planet. The white dwarf is 
ejected and the main-sequence star and planet remain in orbit around the neutron star. As the 
main-sequence star evolves, mass transfer onto the neutron star leads to the formation of the 
currently observed millisecond pulsar. In this case, the orbit of the binary pulsar should have been 
circularized very efficiently by tidal interactions during the mass-transfer phase, leaving a residual 
eccentricity e\ ~ 10~ 4 (Phinney 1992). According to the above results, the planet cannot later 
induce an eccentricity nearly as large as observed today. 

However, Sigurdsson (1993) also points out that because of interactions with passing stars, the 
standard secular perturbation result (eq. [1]) may not apply over a timescale comparable to the age 
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of the system. If each interaction arbitrarily "resets" the initial conditions for the secular pertur- 
bation problem, then over a timescale t p >• 2it / g the eccentricity e\ may grow in a random walk 

1/2 

manner, with e\(t) ~ (Aei) max iV i ^ t , where Ni nt = t/ti nt is the average number of interactions in a 
time t. The mean time between interactions for an object of size 0,2 near the center of M4 (density 
p = 10 4 /O4 Mq pc -3 and velocity dispersion a = 5a$ km s _1 ) is t{ nt ~ 10 8 yr p^ "^(c^/lO au) _1 . This 
is certainly not much shorter than the age of the system, t p <J 10 9 yr (TAT). Even with optimistic 
values for all parameters, we get Ni nt <J 100 which is not sufficient to produce ei(t v ) = 0.025. More- 
over, since the outer orbit for a planet is very "soft" (binding energy much smaller than the typical 
energy of a cluster star), each interaction has a significant probability of disrupting the system. 
Thus a number of interactions as large as 100 would leave a very low probability of survival for the 
triple configuration (TAT also use this likely small value of ti nt /t p to argue against a planet). 



3. Short-Term Perturbations 



Currently available timing data for PSR B1620-26 have been obtained over a time t Q b s ~ 5yr. 
Clearly, we have Pi <C t bs <C P2 for most allowed configurations. Therefore, to a good approx- 
imation, the secular perturbations observable today or over the next few years can be calculated 
assuming that the second companion has a fixed position in space. The components of the perturb- 
ing (tidal) force per unit mass acting on the pulsar are then 

R = Gm 3 2r [3 sin 2 6> 2 cos 2 (</> 2 - v) - 1] (2) 

in the radial direction (r), 

B = j— sin 2 62 cos(</>2 — v) sin(</>2 — v) (3) 



T 2 



in the perpendicular direction (</>), and 



N = j— 3 sin 62 cos 62 cos(</>2 — v) (4) 



r 2 



in the vertical direction (z). Here r2, #2 an d 4>2 are the fixed spherical polar coordinates of the 
second companion (with the origin at the center of mass of the binary pulsar and </>2 measured from 
pericenter in the orbital plane), and r and v are the pulsar's radius and true anomaly. We have 
assumed that r <C ^2- 

Using Lagrange's planetary equations (e.g., Danby 1988), one can calculate the changes in oj\ 
and ei, integrated over one period of the inner orbit. For small e\ we find 

Awi = 3tt 77 [sin 2 6> 2 (5 cos 2 </> 2 - 1) - l] , (5) 

and 

— 15-7T 

Aei = — - — 77 ei sin 2 6 2 sin(2</> 2 ), (6) 
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where rj = [m 2 /(m p + m,i)](ai/r 2 ) 3 . These results are illustrated in figures 3 and 4, where we also 
demonstrate very good agreement with numerical solutions. As an additional check, one can easily 
verify that integrating equation (5) over </> 2 yields the familiar result for the precession of an orbit 
due to a distant circular ring. For a noncoplanar system the perturbation of the inclination is 

3-7T 

Aii = — V sin(20 2 ) cos(wi + </> 2 ), (7) 

which may be detectable through a change in the pulsar's projected semi-major axis a p s'mi. The 
longitude of the node is also perturbed, but this quantity is essentially unmeasurable for a binary 
pulsar. There is no perturbation of the semi-major axis (or orbital period). 

The measured values of the various derivatives of the pulse frequency / can be used to constrain 
the amplitude of the changes predicted by equations (5)-(7). A complete analysis of the problem 
will be presented elsewhere (Rasio 1994). Here, for simplicity, consider the case of a nearly- circular 
outer orbit (e 2 ~ 0) and assume that sinii = sin # 2 = 1. Since f 2 /f <C / and |////| <C |/|, we can 
write the acceleration-induced frequency derivatives as 

• a • n f Gm 2 . 

f = f = 2~ sm(wi + </> 2 ) (8) 

C C C&2 

■■ a • n / G 3 / 2 m 2 Nl / 2 

/ = / = ( m P + mi + m2 > cos(wi + 2 ) (9) 

c c a 2 ' 

■ a • n / G 2 m 2 

/ = /— — = ~ a 5 + m l + m 2) 5in(Ml + (p2) (10) 

Here oj\ = 117 deg is the longitude of pericenter measured from the ascending node (7r/2 — oj\ — </> 2 is 
the angle between the line of sight and the direction to the second companion), a is the acceleration 
of the binary pulsar, and n is a unit vector in the direction of the line of sight. 

For m 2 <C m p + mi, it is straightforward to solve the system of equations (8)-(10) for the 
three unknowns m 2 , a 2 , and </> 2 , given measured values of the three frequency derivatives. In 
general, a few iterations over m 2 may be needed. If we assume that the present value of / is 
determined predominantly by the acceleration of the pulsar (TAT), and adopt the preliminary 
value of / = 10~ 32 s~ 4 reported by Backer (1994), we find the solution m 1 = 1.1 X 10~ 2 M @ , a 2 = 
5.2 X 10 14 cm, and </> 2 = —94 deg. The corresponding rates of change given by equations (5) and (6) 
are tbi = Alj 1 /P 1 = -1.4 X 10~ 4 degyr -1 , e x = Ae 1 /P 1 = -1.1 X 10~ 8 yr -1 . These values should be 
compared to the current (1<t) upper limits derived from the timing data, \u>i \ < 3x 10~ 4 deg yr -1 and 
|ei| < 10 _6 yr _1 (S. Thorsett, private communication). Michel (1994) has calculated more general 
zero-inclination solutions for e 2 7^ 0. In this case the longitude of pericenter w 2 of the outer orbit 
appears as an additional free parameter. For all solutions given by Michel (varying both e 2 and w 2 ), 
we find that equations (5) and (6) predict —1.6 X 10~ 4 degyr -1 < uji < +0.9 X 10~ 4 degyr -1 and 
-0.9 X 10" 7 yr _1 < ei < +0.7 X 10~ 7 yr _1 . The solutions have m 2 varying in the range 10~ 3 M @ to 
1 M & and P 2 between 50 yr and 1300 yr. If we except special values of the angles, we find that most 
solutions give —1.6 X 10~ 4 degyr -1 < < —1.1 X 10~ 4 degyr -1 , which is just below the present 
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upper limit. Thus we expect u>i, in particular, to become measurable very soon. For comparison, 
the general relativistic precession rate is u>\gr = 4.5 X 10~ 5 degyr _1 [(m p + mi)/1.7 M ] 2 / 3 . 

I am very grateful to S. Thorsett for communicating the results of observations in progress 
and for many useful discussions, and to P. Nicholson for a thorough reading of the manuscript. I 
also thank D. Backer, J. Bahcall, C. Bailyn, P. Hut, R. Malhotra, C. Michel, and S. Sigurdsson for 
helpful discussions and comments. The hospitality of the Aspen Center for Physics is also gratefully 
acknowledged. This work has been supported by a Hubble Fellowship, funded by NASA through 
Grant HF-1037.01-92A from the Space Telescope Science Institute, which is operated by AURA, 
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Fig. 1. — Long-term secular evolution of the binary pulsar's eccentricity for two representative cases. 
In (a), the second companion is a planet of mass rri2 = 80 M@ in an orbit of period P2 = 10 yr 
and eccentricity e2 = 0.5. In (b) it is a star of mass rri2 = 0.8 M & , with P2 = 120 yr and e2 = 0.8. 
Coplanar orbits have been assumed, and e\ = at t = 0. In (a), the eccentricity never grows to a 
value comparable to what is observed today (ei = 0.025). 
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Fig. 2. — Maximum eccentricity induced in the binary pulsar's orbit as a function of the eccentricity 
of the second companion's orbit. Conventions are as in Fig. 1, except that e2 is varied. The dots 
are from numerical integrations of the three-body problem, while the solid straight lines show the 
results of analytic perturbation theory. The vertical dashed lines show the dynamical stability 
limits. Note that for the planetary case (a), the maximum induced eccentricity remains much 
smaller than observed today, even when e2 is very close to the limit for stability. 



Fig. 3. — Short-term variations of the longitude of pericenter oj\ and eccentricity e\ for a config- 
uration with m2 = 0.8 M & , ai/r2 = 10~ 2 , </>2 = 0.17T, and sin 62 = 1. See text for details. The 
solid lines show the results of a numerical integration (performed with rri2 on a circular orbit). The 
dashed lines show the analytic results (eqs. [5] and [6], which assume a fixed position for 7722). 



Fig. 4. — Amplitude of short-term secular variations as a function of the longitude </>2 of the second 
companion. All other parameters have the same values as in Fig. 3. The dots are from numerical 
integrations. The solid lines show the analytic results. 



